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The Fredholm equations for one-dimensional two-component Fermions with repulsive and with 
attractive delta-function interactions are solved by an asymptotic expansion for A) strong repulsion, 
B) weak repulsion, C) weak attraction and D) strong attraction. Consequently, we obtain the first 
few terms of the expansion of ground state energy for the Fermi gas with polarization for these 
regimes. We also prove that the two sets of the Fredhom equations for weakly repulsive and attractive 
interactions are identical as long as the integration boundaries match each other between the two 
sides. Thus the asymptotic expansions of the energies of the repulsive and attractive Fermions are 
identical to all orders in this region. The identity of the asymptotic expansions may not mean that 
the energy analytically connects. 
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I. INTRODUCTION 

One-dimensional (ID) Fermi gases with delta- function 
interaction are important exactly solvable quantum 
many-body systems and have had tremendous impact 
in quantum statistical mechanics. The two-component 
delta-function interaction Fermi gas with arbitrary po- 
larization was exactly solved by Yang l] using the Bethe 
ansatz hypothesis in 1967. Sutherland 0] generalized 
the ansatz to solve the ID multicomponent Fermi gas 
with delta-function interaction in 1968. The study of 
multicomponent attractive Fermi gases was initiated by 
Yang _3J and by Takahashi Q in 1970. Since then exactly 
solvable models have been extensively studied by a vari- 
ety of methods developed in the context of mathematical 
physics, see 0, H, @. In particular, recent breakthrough 
experiments on trapped fermionic atoms confined to one 
dimension Q have provided a better understanding of 
significant quantum statistic effects and novel pairing na- 
ture in quantum many-body systems. The observed re- 
sult is seen to be in good agreement with the result ob- 
tained using the analysis of exactly solved models [8l-[T2| . 

Despite the Bethe ansatz equations for the ID two- 
component delta-function interaction Fermi gas with ar- 
bitrary polarization were found long ago [l|, it was not 
until much later that this model began to receive more 
attention from cold atoms jT3j . The asymptotic ground 
state energy of the Fermi gas with polarization was stud- 
ied by the discrete Bethe ansatz equations in a strongly 
and weakly interacting regimes in 1J|. But it turns 



out that the asymptotic expansion of the discrete Bethe 
ansatz equations can only be controllable up to the lead- 
ing order correction to the interaction energy. However, 
the fundamental physics of integrable models are usu- 
ally determined by the set of the generalised Fredholm 
integral equations in the thermodynamic limit, see an 



insightful article by Yang [l5| . The solutions of the Fred- 
holm equations have not been thorou ghly investigated 
analytically except few limiting cases [isl |T6j and nu- 
merical result @,M[i3, EE]- It is usual a difficult task to 
solve those Fredeholm equations analytically. It is highly 
desirable to find a systematic way to treat the generalized 
Fredholm equations. 

In the present paper, we develop a systematic method 
to solve asymptotically the Fredholm equations for ID 
two-component Fermi gas with delta function interac- 
tion and with polarization in four regimes: A) strongly 
repulsive regime; B) weakly repulsive regime; C) weakly 
attractive regime and D) strongly attractive regime. The 
first few terms of the expansion of the ground state en- 
ergy for the Fermi gas with polarization are obtained 
explicitly for these regimes. We also address the analyt- 
ical behaviour of the ground state energy at vanishing 
interaction strength. This method which we develop can 
be directly applied to ID multicomponent Fermi gases. 
The result of ID K-component fermions will be reported 
in the second paper of this study [l^ . 



II. THE FREDHOLM EQUATIONS 



The Hamiltonian for the ID iV-body problem [1|, [2 



H = — 



2m 



N 

E 

1=1 



giD ^ S{xi 

l<i<j<N 



X,). (1) 



describes N Fermions of the same mass m with two 
internal spin states confined to a ID system of length 
L interacting via a (5-function potential. For an irre- 
ducible representation [2^■^, l^t-^^]^ the Young tableau 
has two columns. Where iV-|- and are the numbers 
of Fermions in the two hyperfine levels | t) and | |) such 



2 



that N = N^ + N^. The couphng constant guj can be ex- 
pressed in terms of the mteraction strength c = —2/aiu 
as giD = fi^c/m where qid is the effective ID scattering 



length [21 



Let 2m = H — 1 for our convenience. We de- 
fine a dimensionless interaction strength 7 — c/n for the 
physical analysis, with the linear density n = N/L. For 
repulsive Fermions, c > and for attractive Fermions, 
c < 0. 

The energy eigenspectrum is given in terms of the 
quasimomenta {ki} of the fermions via E = ^ ^j=i ^'j, 
which in terms of the function eh{x) ~ {x + ibc/2)/{x — 
ihc/2) satisfy the BA equations 

exp{ikiL) = Y[ ei {h - A^) , 

AT Afj. 

Yl ei (A„ - kj) = - li 62 (A„ - A^) , (2) 

j=l 13=1 

where i = 1, 2, . . . , TV and a = 1, 2, ... , N^. The param- 
eters {Aq.} are the rapidities for the internal hyperfine 
spin degrees of freedom. The fundamental physics of 
the model are determined by the set of transcendental 
equations which can be transformed to the generalised 
Fredholm types of equations in the thermodynamic limit. 
This transformation was found by Yang and Yang in se- 
ries of papers on the study of spin XX Z model in 1966, 
see an insightful article 



A. Repulsive regime 

For repulsive interaction, it is shown from ([2]) that the 
Bethe ansatz quasimomenta {ki\ are real, but all {Aq} 
are real only for the ground state, see Ref. 0. There 
are complex roots of Aq called spin strings for excited 
states. In the thermodynamic limit, i.e., L,N ^ 00, 
N/L is finite, the roots of the Bethe ansatz equations 
are dense enough in the parameter space. Therefore 
we can define the particle quasimomentum distribution 
function ri{ki) — l/[L{ki~ki^i)] in the quasimomentum 
space. Here ki and kij^i are two conjunction quaisimo- 
mentua. Similarly, the distribution function of the spin 
rapidity is defined r2{Xi) = l/[L{Xi — Aj+i)] in spin pa- 
rameter space. In order to unify the notations in the 
Fredholm equations, we replace the parameter A by fc for 
the distribution function of the spin rapidity. Thus the 
above Bethe ansatz equations ^ can be written as the 
generalized Fredholm equations 



ri{k) 
r2{k) 



1 

2^ 



Kx{k ~ k')r2{k')dk' 



Ki(k - k')ri(k')dk 



K2{k~k')r2{k')dk' 



(3) 



(4) 



The associated integration boundaries Bi , B2 are deter- 
mined by the relations 



n : = N/L = [ ri{k)dk, 
; NJL = I ' r2{k)dk, 



(5) 



-B2 



where n denotes the linear density while is the density 
of spin-down Fermions. The boundary Bi characterizes 
the Fermi point in the quasimomentum space whereas the 
boundary B2 characterizes the spin rapidity distribution 
interval with respect to the polarization. They can be 
obtained by solving the equations in [5] In the above 
equations, we denote the kernel function as 



Ki{x) 



1 



27r (£c/2)2 + a;2 ' 
The ground state energy per unit length is given by 



E = 



k'ri{k)dk. 



(6) 



(7) 



The magnetization per length is defined by Sz = {n — 
2n{)/2. Through the boundary conditions ([S|), the 
ground state energy ([7]) can be expressed as function of 
total particle density n and magnetization Sz- In the 
grand canonical ensemble, we can also get the magnetic 
field h and chemical potential via 



h 



^ dEjn, Sz 
dsz 



dE{n, Sz) 
dn 



(8) 



B. Attractive regime 

For attractive regime, i.e. c < 0, it is found from ([2]) 
that complex string solutions of ki also satisfy the Bethe 
ansatz equations. Thus the quasimomenta {ki} of the 
fermions with different spins form two-body bound states 



k'i ±i^c, accompanied by the real spin 



i.e. 

parameter k^ . Here i — 1 , . . . , N^ . The excess fermions 
have real quasimomenta {kj} with j — 1, . . . ,N ~ 2Ni. 
Thus the Bethe ansatz equations are transformed into 
the Fredholm equations regarding to the densities of the 
pairs p2 (fc) and density of single Fermi atoms pi (k) . They 
satisfy the following Fredholm equations [3|, |^] 

pi(fc) = 1^+ K,{k-k')p2ik')dk' 



2tt 



-Q2 



P2[k) ^ ^ + I Ki{k-k')pi{k')dk' 



+ / K2{k~k')p2{k')dk' 



-B2 



(10) 



Here c < in the kernel functions Ki{x). Here the in- 
tegration boundaries Qi and Q2 are the Fermi points of 



3 



the single particles and pairs, respectively. They are de- 
termined by 



N 
1 



= 2 / p2{k)dk + / pi{k)dk 



Qi 



=: — = p2{k)dk. 
^ J-Q2 

The ground state energy per length is given by 



(11) 



E= {2k^ p2{k)dk+ / k^pi(k)dk. (12) 

J-'Q2 J-Qi 

In a similar way, the magnetic field and chemical poten- 
tial can be determined from the relations ([5]) . In next sec- 
tion, we will discuss solutions and analytical behaviour 
of the Fredholm equations. 



III. ASYMPTOTIC SOLUTIONS OF THE 
FREDHOLM EQUATIONS 



After some algebra, we obtain 

ro(fc) = ^ - |Jc(3), Y2{k) « lc(3). 

Here C(z) is the Riemann's zeta function. Then we obtain 



1 , nln2 'inm i^^2^tl\+o{c-% (16) 



27r TTC 



We see clearly that for strong repulsion the distribution 
of ri(fc) is very flat and it is a constant up to the order 
of 1/c^ correction. This naturally suggests that the ID 
Fermions with strong repulsion can be treated as an ideal 
particles with fractional statistics. Substituting (IT6)) into 
linear density ([S]) and energy ([T]), we obtain 



E 



Bi ( 2nln2 
— 1 + 

TT \ C 

B\ ( 2nln2 



2c3 

3£;c(3) 

2c3 



2c3 
9n37r2C(3) 
10c3 



A. Strong repulsion that give 

The strong coupling condition cL/N ^ 1 naturally Bi = 

gives the condition c ^ Bi, where the Fermi boundary 
Bi oc mr according to the Fermi statistics. For balanced 
case, the numbers of spin-up and -down Fermions are 
equal. Thus there is no finite "Fermi" points in spin 
parameter space, i.e. the boundary B2 — t- 00. Taking a E ~ 

Taylor expansion with the kernel function Ki{k — k') in 
at fc' = 0, we obtained with an accuracy up to the 
order of 1 /c^ 



^ 21n2 / 21n2^^ 8(ln2)3 



^\(3) 



7 



+ 0(c-^), 



41n2 12(ln2)2 32(ln2)= 



(17) 



(18) 



r2(/c) « nKi{k) 



E 
2^ 



3c 



K2ik - k')r2{k')dk' . 



(13) 



Here E is the ground state energy per length. By tak- 
ing Fourier transformation with (|13p. we may obtain the 
distribution function 



r2{uj)= [n — j / [2 cosh — 



(14) 



Substituting (ITil) into the Fredholm equation we 
have 



ri(fc) 



1 1 

2^ ^ 2^ 



e-^"l"lf2(a;)e'"'=da; 



(15) 



where 



r„(fc) 



°° e''^''Lu°'duj 

1 _|_ p\i^\c 



We see that the energy is given in terms of the dimen- 
sionless strength 7 = c/n. The leading order of I/7 
correction was found in [l^, [l^ . Actually, the two sets 
of the Fredholm equations can be converted into dimen- 
sionless units. Therefore the ground state energy can 
be written as analytical functions of 7 except at 7 = 0. 
This ground state energy is a good approximation for the 
balanced Fermi gas with a strongly repulsive interaction 
(good agreement is seen for cL/N > 8) , see Figure [1] and 
Figure [D In these figures, solid lines are obtained from 
the ground state energy ([7]) and (fT2|) with the numeri- 
cal solutions to the two sets of the Fredholm equations 
^ and (HI) for repulsive regime and © and ([TU|) for at- 
tractive regime. The dashed lines are plotted from the 
asymptotic ground state energy for the four regimes. 

However, it seems to be extremely hard to obtain a 
close form of the ground state energy of the gas with an 
arbitrary spin population imbalance in repulsive regime. 
This is mainly because that the distribution function 
r2{k) spans in the region —B2 < k < B2, where the 
integration boundary B2 can vary from zero to infinity 
as the polarization changes. The integration boundary 
B2 decreases as the polarization increases. An intuitive 
way of understanding this point is that zero polarization 



4 



corresponds to ^2 = oo while fully-polarized case corre- 
sponds to B2 = 0. From dressed energy formalism [loj . 
we can easily see this monotonic relation between the 
Fermi boundary and polarization by analysing the band 
filling under external field. For high polarization and 
strong repulsion (i.e., <C N), we have the conditions 
c ^ Bi, B2, that allows us to do the following Taylor 
expansion: 



r2{k) - — 



1 cri{k') 
2^ 



1 4 



1 - 



-2kk' 



fc2 



dk' 



B2 



K2{k - k')r2{k')dk' 



B2 



n 1 



Ki{k)-n^K2{k) + 0{c~^). (19) 



Here we denote = N^/L. We notice that the leading 
order of the distribution function r2{k) is proportional 
to 1/c. Furthermore, taking Taylor expansion in we 
obtain 



1 1 

ri(fc) ~ ^ + ^ 



1 

2^ 



B2 



cr2{k') 



32 4 



1 



^2kk' 



+ 0{c-'). 



fc2 



dk' 
(20) 



From the asymptotic distribution functions (|19p and ([^U)) , 
we calculate the density 
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FIG. 1: (Color online) The ground state energy per length 
vs logarithmic 7 — cL/N in the unit of 2m l'^ • compar- 
ison between the asymptotic solutions and numerical so- 
lutions of the Fredholm equations for polarization P = 
0, 0.2, 0.4, 0.6, 0.8, 1.0. In attractive regime, the biding en- 
ergy Eb ~ — was subtracted. The crossing of the two low- 
est curves in attractive regimes indicates a relation between 
the critical polarization and interaction, where the chemical 
potential of single Fermions excess the chemical potential of 
the pairs. An excellent agreement between our asymptotic 
results and numerical plots is seen for A) strong repulsion, B) 
weak repulsion, C) weak attraction and D) strong attraction 



10' 



ri[k)dk K. — 1 -I- — !^ 



B. Weak repulsion 



that gives 



Bi 



4ni 



16^2 



3? 



(21) 

From the energy ([7]) and the distribution function ri(fc) 
(1201) , we may obtain an asymptotic ground energy of the 
highly polarized Fermi gas with a strong repulsion (c ^ 
Bi, B2) 



E 



1 



c 



487.2 



(22) 



In fact, for strong repulsion, the interacting energy in 
the ground state energy of the highly polarized Fermi gas 
solely depends on the BA quantum number iV^, A struc- 
ture can be found for ID K-component fermions [l9| . By 
numerical checking, we see that for 7 > 8 and polariza- 
tion P ={N^- N^)/{N^ + N^) > 0.5, the energy ^ is 
very accurate, see Fig. O 



For weakly repulsive regime, it is convenient to rewrite 
the Fredholm equations ^ and (|3]) as 



ri(fc) = ^+/ K,{k-k')r2{k')dk', (23) 



r2{k) 



1 

2^ 



Ki{k - k')ri{k')dk' (24) 



|fc'|>Si 



The derivation of is straightforward by using Fourier 
transform of the Fredholm equations ([3]) and (HJ, where 
for our convenience in the study, we actually used 



t(fc) 



r„i(fc) |fc| < B„ 

|fc| > B„ 

(fc) |fc| > B„ 

Ifcl < B^ 



(25) 



with ?Ti = 1, 2 in Fourier transformation. These Fredholm 
equations are valid for arbitrary polarization including 
the balanced case. In the following unification of the 
ground state energy, we assume Bi > B2 as an ansatz. 
In the light of Takahashi's unification of the ground state 
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FIG. 2: (Color online) The ground state energy per length 



vs -y — cL/N in the unit of 



comparison between the 



asymptotic solutions and numerical solutions of the Fredholm 
equations for polarization P = 0, 0.2, 0.4, 0.6, 0.8, 1.0. An 
excellent agreement between our asymptotic expansion results 
and numerical plots is seen in (a) for weakly repulsive and 
attractive regimes, (b) for strongly repulsive regime and (c) 
strongly attractive regime. 



energy [2j], we give the ground state energy per length 

(26) 

Ki{k - k')ri{k')dk' 



n2 1 rB2 

E = -^ + — H{k,Bi)dk 



B2 
-By 



fc'|>-Bi 



H{k,Bi)dk, 



where 



H{x,y) 



-)TTgy{x)+yc 



1 , 4(a; - vf + 

-—XC In — -t; T 

2 4(a; + yY + 



1 - G+(2/,a;), 
tan 



1{x~y) 



± tan 



2(a; + y)' 



From ([5]), we find that the integration boundaries B\ and 
B2 satisfy the following conditions 



r2(fc)G+(Bi,fc)dfc, (27) 
ri{k)G-{k,B2)dk, (28) 





_ Bj_ 


1 


T 


TT 


TT 






1 


L 


TT 


TT 



\k\>Bi 



in this weakly repulsive regime. Using the condition 1 
we may obtain the integration boundary B2 (up to an 
order of c contributions). 



B2 



^ c 4{B,+B2)^+c^ 
+ 4^ 4(Si - B,y + c2 

-i^l^tan- ^ 

TT 

, {Bi + B2) 



2{Bi - B2) 



tan 



2(Bi+B2) 



(29) 



B2 = n.TT + f h, ^^-^ + 0{c^). (30) 
Ztt \Bi - U2\ 

The logarithm term in (1^^ is converged as Bi — i?2. 
However, the logarithm term in (I30p becomes divergent 
as Bi = B2. This divergent term in the ground state 
energy can be canceled out. Here we see a subtlety of this 
asymptotic expansion. Similarly, we calculate the Fermi 
momentum Bi by definition ([5]) and the distribution (j23p 



Bi 



c , 4(Bi+S2)2+c2 

UfTT H In — -:r 

^ 47r 4(Si - B2Y + c? 

c 



-i^l^tan- 



TT 2(Bl - B2) 

c 



^(^l±^tan- 



TT 



2(Bi+B2)' 



(31) 



or 



The ground state energy per length in weakly repulsive 
coupling limit can be expressed in terms of the Fermi 
boundaries 



E 



B? B^ 2c 

— + — + -^BiB2 

'dTT 3Tr TT^ 



(33) 



Substituting (|30|) and (|32|) into (pSj) . we obtain the 
ground state energy of the Fermi gas with a weakly re- 
pulsive interaction and with polarization 



E = -n\TT^ + -uItt^ + 2cn^n4, + 0{c^). (34) 
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This leading order correction to the interaction energy 
indicates a mean field effect. By numerical checking, we 
see that the energy (|34| agrees well with numerical result 
in this weak coupling regime, see Fig. [5J 

For balanced case, the integration boundary B2 tends 
to infinity. It is different from above setting where we 
consider Bi > i?2- Fredholm equations ^ and @ (or 
from and (HH)), can be simplified with the help of 
Fourier transformation 



ri(fc) 



1 

TT 



if2(fc-fc')riout(fc')*'- (35) 



fe'|>Si 



However, for \k\ > Bi, we find 

1 r°° 

^iout(fc) ^ R{k - k')run{k')dk' , (36) 

where the function R(k) is given by 
1 n 1 



i?(fc) = 



27r 1 + e-l'^l'^ 
1 / c c3 
I 4fc2 ^ 8fc4 ^ 



We see clearly that the second term in (p6|) gives a contri- 
bution 0(c). Substituting the leading term ri{k) = 1/2tt 
for the region \k\ > Bi into the distribution function 
ri(fc) (j35[) . we obtain the distribution function ri(fc) for 
|fc| < 



ri(A:) 



1 

TT 



1 

2^ 



tan 



tan 



Bi-k 



From the relation n = J^g_^ ri{k)dk, we find 

1 



1 



47rBi 



■log — ^ 



■ tan 



2Bi 



(3.7) 



Then we obtain the balanced ground state energy for 
weakly repulsive interaction 



E 



37r 



Bi 



1 



c , 4Bl + c^ 



■ tan 



2Bi 



« c 
3W1 



(38) 



It is clearly seen that up to the order O(c^) the ground 
state energy of the balanced gas with a weak repulsion is 
converged as c (or say cL/iV ^ by a rescaling in 
the above equations). By substituting Bi into ([55]). we 
find that the logarithmic term is canceled . The ground 
state energy for balanced case is given by 



E = -^n^T^^ + \n^c + O(c^). 



(39) 



We will further discuss the continuity of the energy at a 
vanishing interaction strength next section. 



Weak attraction 



For weakly attractive regime, the Fredholm equations 
© and are rewritten as 



Pi(fc) = ^ 



P2{k) 



1 

2^ 



Ki{k-k')p2{k')dk', (40) 
Ki{k~k')pi{k')dk'. (41) 



|fc'|>Qi 



These Fredholm equations are valid for arbitrary polar- 
ization. It is seen clearly that the Fredholm equations 
and ([M)) for repulsive regime and the Fredholm 
equations (^0]) and (|41l) for attractive regime are iden- 
tical as long as the integration boundaries match each 
other between the two sides. Similarly, for unification 
of the energy of the gas with a weak attraction, we as- 
sume Qi > Q2. In the above equations, Qi and Q2 are 
determined by 





_ Qi 


1 


T 


TT 


TT 




_ Q2 


1 


L 


TT 


TT 



Qi 

-Qi 
f 

\k\>C 



P2ik)G+iQi,k)dk, (42) 
Pi{k)G-{k,Q2)dk (43) 



that indeed match the integration boundaries Bi and B2 
dill), (IMl) in the region Bi > B2 and Qi > Q2- 

Furthermore, substituting (j40|) and (|4T|) into the 
ground state energy ([T2|) . we obtain 



(44) 

H{k,Qi)dk. 



^ ^ % + ^ I Hik,Qi)dk 



Qi 



Ki{k - k')pi{k')dk' 



k'\>C 



From the equations and (|43l) . we find 



Q2 



{Qi 


-Q2) 




TT 


{Qi 


+ Q2) 




TT 




1^1 m 

47r 


{Qi 


-Q2) 




TT 


{Qi 


+ Q2) 


TT 



4(QiH 


^Q2?- 




4(Qi- 


~Q2Y- 


1- c2 


tan^^ 


\c\ 




2(Qi - 


Q2) 


tan""'^ 


\c\ 




2(Qi -f 


Q2) 


4(QH 


hQ2)2- 




4(Qi- 


-Q2Y- 




tan^^ 






2(Qi - 


Q2) 


tan"^ 






2(Qi + g2) 



O(c^), (45) 



0(c2), (46) 



Indeed, by calculating the ground state energy (|44l) with 
the integration boundaries (|45|) and (j46|) for weakly at- 
tractive interaction, we do find a similar form of the 
ground state energy 



E 



n'\TT + o"-!'^ ~ 2|c|n^n; -I- 0{&'). 



(47) 
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We see that the asymptotic ground state energy (p4)) and 
(|47p continuously connect at c = for arbitrary polar- 
ization, see Figures [T] and [21 

For balanced attractive regime, the Fermi boundaries 
Qi = and Q2 is finite, the Fredholm equations ([9]) and 
(or (go]) and (gl])) reduce to 



P2(fc) = - + / K2{k - k')p2{k')dk' . (48) 

By a iteration, the Fermi boundary Q2 is obtained from 
n = 2 J^Q^ P2{k) in a straightforward way 



Q2 



1 



47rQ2 



log 



4g 



1 -1 
- tan ^ 



49) 



that gives a similar form as the Fermi boundary Bi has, 
see (|37l) . After a length algebra and iteration, we obtain 
the ground state energy 



2 

3^' 



1 



3 /|c| , 4Qi + 



4 

--tan- 



8^ 
3Q2 



+ 0(c2). (50) 



It is clearly seen that up to the oder O(c^) the ground 
state energy of the balanced gas with an attractive inter- 
action is converged too as c — 0~. By substituting Q2 
into ((SO)) , we find that the logarithmic term is canceled 
OTit. Thus the energy is given by 



E=^n'n'-ln'\c\ + 0{c') 



(51) 



that continuously connects to the energy p9p at c — > 0. 
But the identity of the asymptotic expansions may not 
mean that the energy analytically connects because of the 
divergence in the small region c ^ iO and the mismatch 
of the Fermi boundaries associated the two sets of the 
Fredholm equations for both sides. Nevertheless, we see 
that under a mapping 

ri{k) < > piik), r2{k)i 5- P2(fc), ci > c, (52) 



the Fredholm equations ^ and ^ with (H?]), ^ for 
repulsive regime and the Fredholm equations (j40p and 
with , for attractive regime are identical for 
Qi > Q2 and Bi > B2- In the above equations c > for 
repulsive interaction regime and c < for attractive in- 
teraction regime are implied. We also see that the ground 
state energy of the gas with a weakly repulsive interaction 
(p6)) and a weakly attractive interaction (|44)) are unified 
under the mapping (j52p . This unification leads to the 
continuity of the energy for this ploarized gas at vanish- 
ing interaction strength, i.e. c — > 0. Thus the asymptotic 
expansions of the energies of the repulsive and attractive 
Fermions with non-zero polarization are identical to all 
orders in the vanishing interaction strength limits as long 
as the conditions Qi > Q2 and Bi > B2 hold. 



The analyticity of the the energy at c = was dis- 
cussed by Takahashi [2^. Takahashi's theorem states 
that a) the energy function f(n^, tt-j,; c) is analytic on the 
real axis of c when n-^ n^; b) f(n-f, n^] c) is analytic on 
the real axis of c except for c = when n-^ = n^. This 
theorem appears not to be true for the region Bi < B2 
and Qi < Q2 in our study. Takahashi's proof of this the- 
orem relies on his Lemma 2, i.e. the function /, density 
n and density of spin-down Fermions ni are analytic as 
functions Q,B and c except for the region c = 0, and 
Q < B. Here Q and B are two integration boundaries. 
Even the identity of the asymptotic expansions of the en- 
ergy may not mean the energy analytically connects due 
to the divergence of the two sets Fredholm equations in 
the limit c — ;> iO and the mismatch of the intervals for the 
density distribution functions. Although we unified the 
two sets of Fredholm equations (HH) and (gO]), (|4T|) 
for arbitrary polarization, the integration boundaries be- 
tween the two regimes are mismatched for the regions 
Qi < Q2 and Bi < B2, i.e. 



T 



TT 



/Si<|fc|< 
B2_l 
TT TT 



r2{k)G+{Bi,k)dk 



r2{k) 



1 



i--G_(fc,Bi; 



A;|>B2 

ri{k) 



ri{k)G-{k,B2)dk 
1 



1 



-G+{B2,k) 



lBi<\k\<B2 

for weakly repulsive regime, and 



dfc(53) 



dfc(54) 



t 



L 



Qi 

TT 



Qi 



P2ik)G+iQi,k)dk 
1 



P2{k) 



'Ql<\k\<Q2 
h 1 



1 



G-(fc,Qi) 



TT 



J\k\>Q2 

Pl(fc) 

Qi<|fc|<Q2 



Pi{k)G-[k,Q2)dk 

1 



1 



:G+(Q2,fc) 



dfc(55) 



dA:(56) 



for weakly attractive regime. It is obvious the signs in the 
last term in each equation are mismatched. In the above 
equations c > for repulsive interaction regime and c < 
for attractive interaction regime are implied. This mis- 
match is clearly seen from the balanced case: -Bi — > Q2, 
B2 ^ 00 and Qi 0. Thus we see that the Fredholm 
equations can not be unified for the region Qi < Q2 and 
Bi < B2 in the vanishing interaction strength. 



D. Strong attraction 

In recent years, strong attractive Fermi gas has been 
received a considerable attention from theory and experi- 
ment due to the existence of a novel pairing state. For the 
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spin-1/2 Fermi gas with strongly attractive interaction, 
two Fermions with different spin states can form a tightly 
bound pair. For the ground state, the model has three 
distinguished quantum phases, i.e. fully paired phase 
with equal number spin-up and -down Fermions, fully po- 
larized phase of single spin-up Fermions and a partially 
polarized phase with both pairs and excess Fermions. 
The key features of this T = phase diagram of the 
strongly attractive spin-1/2 Fermi gas were experimen- 
tally confirmed using finite temperature density profiles 
of trapped fermionic ^Li atoms 

Here we calculate the ground state energy (|12p from 
the Fredholm equations ^ and (|10p with the integration 
boundaries Qi,Q2 that characterize the Fermi points of 
two Fermi seas, i.e. the Fermi seas for excess Fermions 
and pairs, respectively. Therefore, for strong attraction, 
i.e. \c\L/N ^ 1, all integration boundaries are finite, 
i.e. Qi and Q2 are finite. In this regime, the condi- 
tions c 3> (5i,Q2 hold for arbitrary polarization. From 
the following calculation, we will see that the conditions 
Qi > Q2 and Qi < Q2 do not change the expression of 
the energy. Therefore, the following result is valid for ar- 
bitrary polarization, including the balance case. In this 
regime, it is convenient to use a notation |c| instead of a 
negative value of c. The ground state energy is calculated 
in the following way 



-Q2 



-tan-^(^+tan-^(^ 



2^"^"^"4(fc- 



|c| 

l)2 + c21 



?l)2 + c2 



dk. (57) 



Furthermore, we consider strong coupling expansion in 
the energy ((57|) . here we assume |c| ^ Qi,Q2- We collect 
contributions up to the order of l/|c|'^, i.e. 



E 



Of 
Svr 



4ni 



32Ql 



5 c- 



Stt c ■ 



P2{k)k^dk. 



(58) 



with 



° Stt 
2 



4n 



W , 32Qi 



|c| 



5|c|^ 37r|cP 
P2{k)k'^dk, Eb 



-Q2 



y (60) 



We calculate Qi from (fTTI) with the density 




Ki{k - k')p2{k') dk 



An^ , IQQin^ , 32Q^ 



|c| 3|c|3 37r|c|3 
Then we obtain the Fermi momentum 




|c| 



I2nl 



(61) 



Similarly, we calculate Q2 from pT|) with the distribu- 
tions ® and (Uni) 



■,2^2/ 



2?it - (2n^ - n^)^ (2n^ - nj,) 



|c| 



+ 



12|d3 



,3 + 32(„^ _ „^)3 



12 c 



62) 



We observe that the kernels in the Fredholm equations 
© and ([TUl) converge quickly with the distribution func- 
tions as the interacting strength |c| increases. This allows 
one to take a proper Taylor series expansion in the dis- 
tribution functions. In this way, from Eq. (|10p we may 
obtain 



P2{k) 
1 

~2^ 



1 



c2 + fc2 2(c2 + A:2)2 

\c\{nt-n^) , |c| 



f + 



2g 



Qi 



pi{k')k'^dk' 



TT 37r2|c|'^ 37r2|cP TT c2 -|- 



27r 



fc2' 



(63) 



In the last equation of the first part in the square 
bracket is the kinetic energy of excess single atoms in- 
cluding marginal interference effect between the single 
atoms and molecules of two-atoms. The second term is 
the total binding energy of the bound pairs. The last 
term characterizes the total energy of the molecules of 
two-atom. We now calculate the Fermi momenta Qi and 
Q2 and the energy of the molecules of two-atom. For our 
convenience, we denote 



E = E^ + E'o + n^Sb 



(59) 



Substituting §^ into the energy E^ (pH) , 

4Q| 8g|_ 32QfQ| 
3tt 97r2|c|3 97r2|c|3 
2ni f'^^ \c\k'^dk (n^ - n^) 



4Q| 
3tt 



,^ + k"^ IT 



37r|cP 



3(8nt - 7n^)Ql 
5|c|3 



■2 4 



fc2 



(64) 
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Substituting equations ([6T|) and (|62l) into the ground 
state energy ([55]) and we obtain the ground state 

energy of the gas with a strongly attractive interaction 
and with an arbitrary polarization 



8n,, 



48n2 



4807if 



5ny) 



(65) 



6 



1 



15|c| 

+15n|^ - UQn^ + 63Tr'^nln-f - eOTr^Ti^Ti^)] . (66) 

We define polarization P = {N^ — Ni)/N = (n-f- — n^)/n, 
then the energy in terms of polarization is given by 



E 



(f-P)72 7r2(l-3P + 3P2 



f5P3 



2m [ 4 
7r2(l-P)(l + P 



48 



67P^^ 



48|7| 

7^2(1 _ p)2(i + 5P + 3p2 ^_ 247P3) 



5472 



7r2(l-P) 



[-15 + 31125P'' + 186l7r2p^ 



1440|7|3 

-15765P^ - 6597r2p4 + 3467r2p3 - Mtt'^P'^ 
WP + tt^ - 105P - 150P2 - 15090P3] } . (67) 

that agrees with the result derived from dressed energy 
equations [l^, [HI ■ This result is highly accurate as be- 
ing seen in Fig. [1] and Fig. [21 From the energies 
([55)1 and ((M)) . we see the bound pairs have tails and 
the interfere with each other. But, it is impossible to 
separate the intermolecular forces from the interference 
between molecules and single Fermions. If we consider 
ni ^ X = — n^^ the single atoms are repelled by the 
molecules, i.e.. 



Where 



^ + 6"^ 



Ax 



122:(a; + n^) 



6 



nV 1 



2% 
Id 



3n2 



In addition, the phase diagram and magnetism can be 
work out directly from the relations ([5]) with the ground 
state energy for the four regimes. The phase boundaries 
of the full phase diagrams may be analytically and nu- 
merically obtained by imposing the conditions — 0, 0.5 
in the conditions ([5]), which has been discussed in litera- 
ture [M3,[i3- 



IV. CONCLUSION 

In conclusion, we have presented a systematic method 
to derived the first few of terms of the asymptotic expan- 
sion of the Fredholm equations for the spin- 1/2 Fermi 
gas with repulsive delta-function and attractive delta- 
function interactions in the four regimes: A) strongly 
repulsive regime; B) weakly repulsive regime; C) weakly 
attractive regime and D) strongly attractive regime. We 
have obtained explicitly the ground state energy of the 
Fermi gas with polarization in these regimes, see the key 
resuh (UHl), (gH), (ED), (gZl), dlS]) and (Uni). By numer- 
ical checking, these asymptotic ground state energy are 
seen to be highly accurate in the four regimes. In weakly 
attractive and repulsive regimes, the ground state ener- 
gies, integration boundary relations and the associated 
two sets of the Fredholm equations have been unified. 
The two sets of the Fredholm equations can be identical 
as long as the associated integration boundaries match 
each other between the two sides. This suggests that 
the asymptotic expansions of the energies of the repul- 
sive and attractive Fermions are identical to all orders 
in this region as c 0. The identity of the asymptotic 
expansions may not mean that the energy analytically 
connects due to the divergence of the two sets Fredholm 
equations in the limit c — iO and the mismatch of the 
associated integration boundaries between two sides at 
some intervals, e.g.. Pi < P2 and Qi < Q2- 

Moreover, the explicit result of the ground state energy 
obtained provides facilities to study universal nature of 
many-body phenomena. The local pair correlation for 
opposite spins can be calculated directly from the ground 
state energy by 



42(0) 



1 



dE{n, Sz)/dc. 



This naturally gives the ID analog of the Tan's adiabatic 
theorem [23| through the relation 



C 



''ID 



where C is called the universal contact, measuring the 
probability that two fermions with opposite spin stay to- 
gether. It was shown [^J] that the momentum distri- 
bution exhibits universal C/k'^ decay as the momentum 
tends to infinity. The significant feature of Tan's uni- 
versal contact is that it can be applied to any many- 
body system of interacting bosons and fermions in ID, 
2D and 3D [2^. In addition, the explicit forms of 
the ground state energies in the four regimes can be used 
to determine magnetism and phase diagram of the sys- 
tem in the grand canonical ensemble. It can help con- 
ceive quantum statistical effect by a comparison between 
the ground state energies of ID delta-function interact- 
ing fermions and spinless bosons. These provide a precise 
understanding of many-body correlations and quantum 
magnetism in the context of cold atoms. The method 
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which we have developed in this paper can be generahzed 
to study ground state properties of ID multicomponent 
Fermi and Bose gases with delta-function interaction. We 
will consider it in 1^. 
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